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ABSTRACT 

The  exact  solution,  the  Bom  approximation,  and  its  variational  improvement  are 
obtained  for  the  scattering  of  electromagnetic  waves  from  a  random  ensemble  of 
systems,  each  consisting  of  two  Rayleigh  cylinders.  The  cylinders  are  parallel,  of  in¬ 
finite  length,  and  of  equal  radius.  Their  separation  varies  randomly  among  ensemble 
members  except  that  the  cylinders  cannot  overlap.  The  intent  is  to  test  a  recently  de¬ 
veloped  vector  stochastic  variational  principle.  The  exact  solutions  are  obtained  for 
the  average  differential  scattering  cross  sections  of  both  the  transverse  electric  (TE) 
and  transverse  magnetic  (TM)  fields  relative  to  the  cylinder  axes  with  normal  plane 
wave  incidence.  The  corresponding  variational  approximations  are  obtained  using  a 
recently  reported  computational  alternative  to  the  more  familiar  dyadic  Green’s 
function  solution.  They  are  in  essential  agreement  with  the  exact  TE  and  TM  solu¬ 
tions,  whereas  the  Born  results  arc  not.  In  particular,  the  variational  results  accur¬ 
ately  account  for  multiple  scattering,  which  is  significant  in  the  exact  TE,  but  not 
TM,  solution,  and  also  account  for  the  difference  in  geometric  polarizability  be¬ 
tween  the  two  solutions. 
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ILLUSTRATIONS 

1  The  scattering  problem  consists  of  two  nonoveriappling,  parallel,  dielec¬ 
tric  cylinders  illuminated  by  a  plane  wave  normal  to  their  axes  and  form¬ 


ing  a  fixed  angle  0  with  the  plane  in  which  their  axes  lie .  9 

r  2  Illustration  of  parameters  that  relate  to  /?,  using  the  cosine  law .  16 
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1.0  INTRODUCTION 


The  potential  advantages  of  variational  approx¬ 
imations  in  scattering  problems  are  well  known.14 
First,  by  virtue  of  variational  invariance,  errors  in  the 
initial  trial  approximation  for  the  fields  on  or  within 
the  scatterer  do  not  lead  to  corresponding  first-order 
errors  in  the  calculated  far  field  scattering  amplitude 
T.  Thus,  good  initial  approximations  can  lead  to  im¬ 
proved  results,  and  estimates  of  the  accuracy  of  vari¬ 
ational  results  can  be  made  in  some  cases.4  4  Because 
the  usual  Schwinger  variational  form  for  the  scatter¬ 
ing  amplitude  is  a  ratio  of  integrals,-  difficulties  in 
computing  the  statistical  moments  of  this  ratio 
resulted  in  very  limited  application  of  this  method  to 
random  scattering  problems  until  the  mid-seventies. 
In  1977,  Hart  and  Farrell2 * 4  demonstrated  that,  for  ar¬ 
bitrary  scatterer  statistics,  the  integrals  whose  ratio 
comprises  the  deterministic  variational  expression 
can  be  averaged  individually  and  recombined  to  form 
an  invariant  expression  for  the  average  scattering 
amplitude  <f>  .  This  form,  i.e.,  <D  =  </V, ) 
x  </V2  >/<D>,  is  inherently  simpler  to  evaluate  than 
the  direct  average  of  the  deterministic  expression  <  T > 
=  <A/,  N2/£>>,  where  N] ,  N2,  and  D  are  integrals  in¬ 
volving  the  fields  at  the  scatterer. 

Hart  and  Farrell  demonstrated  this  result  for 
<  7>and  <  ITI2 )  in  the  case  of  scalar  wave  scattering 
from  objects  or  surfaces  with  homogeneous  boun¬ 
dary  conditions.  An  extension  was  made  recently5 *'' 
to  vector  wave  scattering  and  inhomogeneous  boun- 


1 L.  Cairo  and  T.  Kahan,  Variational  Techniques  in  Electromag¬ 
netism.  Gordon  and  Breach,  Science  Publishers,  Inc.,  New  York 
(1965). 

2  P.  M.  Morse  and  H.  Feshbach,  Methods  of  Theoretical  Physics, 
Parts  I  and  II,  McGraw-Hill,  New  York(!953). 

JD.  S.  Jones,  “A  Critique  of  the  Variational  Method  in  Scatter¬ 
ing  Problems,”  IRE  Trans.  Antennas  Propag.  AP-4,  297-301 
(I9J6). 

4R.  W.  Hart  and  R'  A.  Farrell,  "A  Variational  Principle  for 
Scattering  from  Rough  Surfaces,”  IEEE  Trans.  Antennas  Pro¬ 
pag.  AP-2S,  708-710(1977). 

5J.  A.  Krill  and  R.  H.  Andreo,  “Vector  Stochastic  Variational 

Principles  for  Electromagnetic  Wave  Scattering,”  IEEE  Trans. 

Antennas  Propag.  AP-28,770-776(l980). 

*R.  H.  Andreo  and  J.  A.  Krill,  “Vector  Stochastic  Variational 

Expressions  for  Scatterers  with  Dielectric,  Conductive,  and 
Magnetic  Properties,"  J.  Opt.  Soc.  Am.  71, 978-982  ( 1981 ). 


dary  conditions.  Specifically,  we  have  derived5  '’  gen¬ 
eral  vector  stochastic  variational  expressions  for  the 
statistical  moments  and  probability  density  functions 
of  arbitrary  polarization  components  of  the  vector 
scattering  amplitude,  e,  ■  T  and  the  differential  scat¬ 
tering  cross  section  le,  ■  Tl:.  These  invariant  ex¬ 
pressions  are  applicable  to  random  scatterers  with  ar¬ 
bitrary  inhomogeneous  and  anisotropic  permittivity, 
conductivity,  and  permeability. 

Of  present  interest  is  the  extent  to  which  the  gener¬ 
al  vector  stochastic  variational  principle,  by  virtue  of 
first-order  error  cancellation,  can  approximate  polar¬ 
ization-dependent  interactive  phenomena,  i.e.,  inter¬ 
ference  and  multiple  scattering,  even  if  the  initial 
field  approximation  does  not.  This  question  is  ad¬ 
dressed  in  this  report  by  calculating  the  exact  solution 
to  a  random  scattering  problem  exhibiting  these  ef¬ 
fects  and  then  comparing  the  result  with  the  corre¬ 
sponding  Born  approximation  and  its  variational  im¬ 
provement. 

Our  approach  is  based  on  a  similar  earlier 
analysis.7 *  *  In  an  effort  to  examine  the  accuracy  and 
efficacy  of  the  original  scalar  stochastic  variational 
principle,4  Gray,  Hart,  and  Farrell’  calculated 
closed-form  variational  and  first-order  perturba- 
tional  approximations  to  transverse  magnetic  (TM) 
wave  scattering  from  a  classic,  random,  perfectly 
conducting  surface9  on  which  Dirichlet  boundary 
conditions  are  satisfied.  The  surface  consisted  of 
many  parallel,  nonoverlapping,  hemicylindrical  pro¬ 
trusions  on  an  infinite  plane.  The  hemicylinders  were 
of  equal  radii,  infinite  length,  and  random  separa¬ 
tion.  The  variational  approximation  for  <  ITI2)  was 
found  to  include  the  sum  of  independent  hemicyl- 
inder  scattering  contributions  as  well  as  a  correction 
term  proportional  to  the  fractional  area  of  the  plane 


7E.  P.  Gray,  R.  W.  Hari,  and  R.  A.  Farrell,  "An  Application  of 
a  Variational  Principle  for  Scattering  by  Random  Rough  Sur¬ 
faces,"  Radio  Set.  13,  333-348  ( 1978). 

8J.  A.  Krill  and  R.  A.  Farrell,  "Comparisons  Between  Varia¬ 
tional,  Perturbational.  and  Exact  Solutions  for  Scattering  from 
a  Random  Rough  Surface  Model,"  J.  Opt.  Soc.  Am.  68. 
768-774(1978). 

9  V.  Twersky,  "Multiple  Scattering  of  Radiation  by  an  Arbitrary 
Planar  Configuration  of  Parallel  Cylinders  and  by  Two  Parallel 
Cylinders,"  J.  Appl.  Phys.  13,  No.  4, 407-414  ( 1952). 
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covered  by  the  hemicylinders.  Comparison  of  this  re¬ 
sult  with  first-order  perturbation  theory  showed  a 
discrepancy  that  persisted  even  when  the  Rayleigh 
limit  was  taken.  Krill  and  Farrell"  were  able  to  obtain 
the  exact  solution  for  the  simplest  multiple  scattering 
case  for  this  perfectly  conducting  surface,  consisting 
of  only  two  randomly  separated,  nonoverlapping, 
Rayleigh  hemicylinders.  They  showed  from  the  exact 
solution  that  the  variational  result  accounted  for 
multiple  scattering  whereas  the  first-order  perturba¬ 
tion  approximation  did  not,  thus  offering  an  ex¬ 
planation  for  the  discrepancy. 

The  present  study  differs  from  the  earlier  one’  *  as 
follows.  First,  we  are  testing  the  recent  extension  of 
the  stochastic  variational  principle  to  (vector)  electro¬ 
magnetic  wave  scattering  from  objects  or  surfaces 
with  generally  nonvanishing  boundary  conditions. 
The  scattering  model  to  be  considered  involves  di¬ 
electric  cylinders.  Because  there  are  differences  be¬ 
tween  the  exact  transverse  electric  (TE)  and  TM  fields 
(including  interactive  effects),  we  will  investigate  the 
extent  to  which  the  vector  variational  principle  ac¬ 
counts  for  these  polarization  phenomena.  The  model 
consists  of  an  ensemble  of  systems,  each  of  which 
contains  two  dielectric,  parallel,  Rayleigh  cylinders 
with  infinite  lengths  and  equal  radii.  A  plane  wave  is 
incident  perpendicular  to  the  cylinder  axes  and  makes 
a  fixed  angle  relative  to  the  plane  in  which  the  axes 
lie.  The  cylinders  are  randomly  separated  but  hard; 
i.e.,  they  cannot  merge  or  overlap,  and  their  refrac¬ 
tive  index  goes  discontinuously  to  the  value  1  at  their 
surface. 


For  normal  incidence  and  arbitrary  incident  polar¬ 
ization,  the  scattered  field  can  be  written  as  the 
superposition  of  the  (decoupled)  TM  and  TE  solu¬ 
tions.""'  The  exact  TE  and  TM  solutions  for  the  dif¬ 
ferential  scattering  cross  section  are  obtained  in 
Chapter  2  for  two  identical  dielectric  Rayleigh 
cylinders  with  axis  separation  greater  than  their  di¬ 
ameter.  In  Chapter  3,  these  solutions  are  ensemble 
averaged  over  a  separation  distance  subject  to  the  no¬ 
overlap  constraint.  Chapter  4  presents  calculations 
using  the  vector  stochastic  variational  principle  to 
improve  the  Born  approximation  for  random  cylin¬ 
der  scattering  and  arbitrary  polarization.  For  this 
simple  random  scattering  problem,  the  average  of  the 
deterministic  expression  could  also  be  calculated  and 
is  presented.  In  Chapter  S,  the  Born  and  variational 
results  are  compared  to  the  exact  TE  and  TM  solu¬ 
tions,  with  particular  attention  to  the  accuracy  with 
which  the  approximations  account  for  polarization 
effects  and  multiple  scattering.  It  is  concluded  that 
the  variational  approximation  accurately  describes 
both  the  polarization  dependence  and  the  multiple 
scattering  found  in  the  Rayleigh  limit  of  the  exact 
solution. 


"'H.  C.  Van  de  Hulsl,  Light  Scattering  by  Small  Panicles,  Wiley, 
New  York,  pp.  297-301, 306-307  (1957). 

nG.  Olaofe,  "Scattering  by  Two  Cylinders,"  Radio  Sci.  5,  No. 
It,  1351-1360(1970). 
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2.0  EXACT  TE  AND  TM  SCATTERING  CROSS  SECTIONS 


Figure  1  presents  the  scattering  configuration  in 
which  two  parallel,  infinitely  long  cylinders,  labeled 
cylinder  + 1  and  - 1 ,  are  illuminated  by  a  plane  wave 
incident  normal  to  their  axes,  which  are  located  at 
i<+,)  and  f'11,  respectively.  The  cylinders  have 
relative  refractive  index  m  and  radius  a,  and  their 
axes  are  along  z.  The  incident  plane  wave  has  electric 
field  E,  =  Ae,  exp(ik,  •  r),  where  k,  is  the  wave  vec¬ 
tor,  A  is  the  amplitude,  and  e,  is  the  polarization  vec¬ 
tor.  The  wave  vector  is  in  the  direction  x,  i.e.,  k,  s 
for,  and  forms  a  fixed  angle  0  relative  to  the  plane  in 
which  the  cylinder  axes  lie.  We  are  interested  in  the 


Ftgur*  1  —  The  scattering  problem  consists  of  two 
nonoverlapping,  parallel  dielectric  cylinders  illuminated  by 
a  plane  wave  normal  to  their  axes  and  forming  a  fixed  angle 
6  with  the  plane  in  which  their  axes  lie.  The  infinite-length 
cylinders’  separation  f  =  f'*1’  -  jJ_1)  is  a  random 
variable,  and  TE  and  TM  scattering  cross  sections  are 
calculated  in  the  Rayleigh  limit. 


differential  scattering  cross  section  I7'l2,  where  T 
ae,  •  T,  for  arbitrary  polarization  e,.  Because  the 
plane  wave  is  incident  normal  to  the  cylinder  axes,  ar¬ 
bitrary  polarization  can  be  treated  as  the  sum  of 
decoupled  TE  and  TM  field  contributions.  The  exact 
TE  and  TM  solutions  will  be  obtained  for  the 
Rayleigh  limit10  mka  and  ka  <  <  1. 

The  exact  solution  for  scattering  from  two  dielec¬ 
tric  cylinders  with  arbitrary,  but  fixed,  axis  separa¬ 
tion  -  f*'11  has  appeared  in  several 

sources.9,1114  As  in  our  two-hemicylinder  solution,8 
we  again  follow  the  procedures  described  by 
Olaofe,1 1,12  which  are  based  on  scalar  wave  functions 
related  to  the  TE  and  TM  fields.10  First,  the  wave 
function  used  to  describe  the  total  field  external 
to  the  scatterers,  is  written  as  the  sum  of  an  incident 
plane  wave  and  scattered  wave  contributions  from 
each  cylinder,  i.e.,  \pex  =  ^inc  +  i/v-  The  incident 
plane  wave  is  expanded  in  a  Bessel  series  in  terms  of 
coordinates  centered  at  they-th  cylinder  axis,  where  j 
=  + 1  or  - 1  (Fig.  1).  The  scattered  field  from  each 
individual  cylinder  is  expressed  as  a  sum  of  cylin- 
dricaily  outgoing  waves  centered  at  that  cylinder  with 
as-yet  undetermined  coefficients.  Finally,  the  addi¬ 
tion  theorem  is  applied  that  expresses  the  outgoing 
waves  from  one  cylinder  m  terms  of  coordinates  cen¬ 
tered  at  the  other  cylinder.  The  resulting  equation  for 
the  external  wave  function,  expressed  in  the  coor¬ 
dinates  centered  on  they-th  cylinder,  is" 

=  A+  fj  em#u)  [/" WjtMJn  ( krU) ) 


+  ,b„H"Hkr^)  +  J„(kr"  )(-!)" 

x  t  rWWl-D1).  0) 


12G.  Olaofe,  “Scattering  Cross  Section  for  Two  Cylinders,”  IEEE 
Trans.  Antennas  Propag.  AP-H,  No.  6, 823-825  (1970). 

,JV.  Twersky,  “Scattering  of  Waves  by  Two  Objects,”  Electro¬ 
magnetic  Waves,  R.  E.  Langer,  ed.,  The  University  of  Wiscon¬ 
sin  Press,  Madison,  pp.  361-389(1962). 

i4N.  Zitron  and  S.  N.  Karp,  “Higher-Order  Approximations  in 
Multiple  Scattering,  1,  Two-Dimensional  Scalar  Case,"  J. 
Math.  Phys.  2,  No.  3, 394-402(1961). 
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where,/'  is  - 1  when  /is  + 1  and  vice  versa;  £  equals  ( 
or  n,  respectively,  when  /  equals  - 1  or  +1;  and  w, 
=  exp  [/Arf01  cos  d] .  The  factor  wy  accounts  for  the 
phase  shift  of  the  incident  wave  that  occurs  when  the 
coordinate  origin  is  translated  to  the  axis  of  cylinder 
/.  The  sum  over  the  first  term  in  the  brackets  is  the 
incident  plane  wave,  the  sum  over  the  second  term  is 
the  scattered  wave  from  the  y-th  cylinder,  and  the 
doubl'  sum  represents  that  from  the/'-th  cylinder. 
Het .  ie  incident  scalar  wave  function  is  assumed  to 
be  a  pivne  wave  of  amplitude  A , .  This  amplitude  can 
be  related10  to  the  corresponding  amplitude  A  of  the 
incident  plane  wave  electric  field.  Because  the  ampli¬ 
tude  of  the  incident  wave  cancels  in  the  final  expres¬ 
sion  for  T  (c.f.  Eqs.  6a,  6b,  40,  and  41),  it  will  not  be 
discussed  further. 

An  expression  for  the  coefficients  in  Eq.  1  can  be 
obtained  by  applying  the  appropriate  boundary  con¬ 
ditions  for  the  fields  at  the  cylinders’  surfaces,10  " 

A  =  -cjfwyew  +  (-l)" 

x  E  j (2) 

( =  —  OD 

n  —  0,  ±1,  ±2, ...,  ±oo, 

where  c„  is  the  n-th  order  scattering  coefficient  for  a 
single  cylinder.10  "  These  latter  coefficients  will  be 
designated  by  c„TE  and  c™  for  the  TE  and  TM  solu¬ 
tions,  respectively. 

For  nonoverlapping  Rayleigh  cylinders  (f  >  2a, 
mka  <  <  1 ,  and  ka  <  <  1),  it  can  be  shown*12  "15  that 
the  infinite  set  of  coupled  equations  in  Eq.  2  may  be 
truncated.  This  conclusion  holds  even  in  the  limit 
that  the  cylinders  touch,  i.e.,  f  —  2a(l  +5,),  where 
6„  —  0.  Series  convergence  considerations  and  ex¬ 
amination  of  c„  reveal*'14'15  that,  in  the  Rayleigh 
limit,  the  A  coefficients  that  are  lowest  order  in  ka 
correspond  to  n  -  0  for  TM  fields  and  n  =  ±  1  for 
TE  fields.  Moreover,  the  series  in  f  may  be  cor¬ 
respondingly  truncated  so  that 


,5Y«.  A.  Ivanov,  Diffraction  of  Electromagnetic  Waves  on  Two 
Bodies,  Nauka  i  Tekhnika  Press,  Minsk,  pp.  66-69  (1968) 
(English  translation). 


A™=-c0™wy,  (3) 

where c0™  =  ix(m2  -  l)(Ara/2)2,  and 
t/cte 

,b\\  =  '  i  [wye*i,J  +  pWj  e*'"  ],  (4) 

where  c,TE  =  -  i*(ka/2)2(m2  - 1  )/(m2  +  1),  and  the 
interaction  parameter  p  =  c^H2(kfi  appears  in  both 
the  exact  and  the  variational  solutions  for  the  TE 
case.  We  note  that  in  the  Rayleigh  limit  of  the  two- 
cylinder  problem,  multiple  scattering  is  expected  to 
be  significant  only  when  the  cylinders  are  in  proximi¬ 
ty.  Observe  in  Eqs.  3  and  4  that  the  only  variable  that 
has  such  a  characteristic  is  p,  which  appears  only  for 
the  TE  case. 

Using  the  standard  relationship10  between  the 
fields  and  the  wave  function,  one  can  show  that  the 
scattering  amplitude  definition. 


results  in  the  following  TE  and  TM  amplitudes: 

f™  =  f  ft™  +  /,TMe-i*rcos(#J-Jl)e-i*,J'l,cos((>,- 

(6a) 

where  it  is  recognized  that 

p*"  =  9  l)  &  es  -  &  m  e, 

r4*"  s  r4  ",  and 
i4’"  _  ft-')  =  -fcosffl,  -  0), 
and 


jte  _  |  B  +  B e~ '4'costw,  “ije  ,4i*  "lomuj-id 

(6b) 

where  tB  =  HjbV\  exp(  -  i(6s  -  (3))  -  yft,TE 
exp(  +  H0,  -  0)),  and  we  have  used  a  coordinate  sys¬ 
tem  centered  on  the  origin  in  Fig.  1 .  The  cross  sec¬ 
tions  1 7™  I2  and  IT56 12  are  readily  obtained  by  tak¬ 
ing  the  absolute  square  of  Eqs.  S  and  6  and  inserting 
Eqs.  3  and  4,  respectively. 
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3.0  ENSEMBLE  AVERAGES  OF  EXACT  SOLUTIONS 


Ensemble  averages  of  the  exact  cross  sections  over 
random  cylinder  separation  f  will  be  evaluated  in  the 
manner  of  Krill  and  Farrell.8  The  cylinder  axes  point 
in  the  z-direction  and  lie  in  the  plane  that  forms  an 
angle  0  relative  to  the  incident  wave.  The  axes  are  lo¬ 
cated  within  the  limits  -  L/2  <  <  +  L/2,  j  - 

±1,  and  are  constrained  not  to  overlap,  i.e.,  f  >  2a. 
The  two-center  probability  density  function  is 


between  the  cylinders  caused  by  their  nonoverlap 
condition. 

The  TE  case  is  treated  by  noting  that  performing 
the  absolute  square  of  Eq.  6b  and  applying  Eq.  4 
leads  to  terms  that  are  of  the  form 


p"p*  cos(yft 

11  -  p2 12  ’ 


(9) 


p\ J< -»]  = 


(TT^iff>2fland 

—  L/2  <+Z./2, 

j  =  ±  1 

0  otherwise 


where  (n, I)  =  0,  1 ,  and  2,  y  =  [A:  cos  0,  k  cos  0,  Ar(cos 
0  ±  cos  9),  or  0],  and  the  interaction  parameter  is, 
again,  p  =  c,  which  approaches  (Ma'/f) 

for  small  Arf  where  M  m  [m2  -  l)/(m:  +  1).  These 
averages  are  evaluated  analytically  in  the  same  man¬ 
ner  as  in  Ref.  8  (see  Eqs.  18  and  20a  of  that  reference) 
to  give 


The  average  <  >  of  an  arbitrary  function  g(  f)  that 
depends  only  on  f  a  -  J"  "  can  be  shown  to 
be8 

<*<»>  =  (L  -2a)2  11  (L  ~  nt*(n  + 

(7a) 


v  v  MU  y2<"  ‘  '  1  * 

<y<f)>  =  (- ir  ^  Va?]  jy— ^  dy.  do) 

where  y  =  \rM  a /f.  The  validity  of  Eq.  10  can  be 
verified  by  demonstrating"  that  it  differs  from  Eq.  9 
by  terms  that  are  higher  order  in  ka.  Equation  10  can 
be  evaluated  analytically16  to  give  (to  first  order  in  i> 
and  second  order  in  M) 


which,  in  the  limit  L  —  oo,  goes  to 8 


<*W>  =i  }jl  -  t/Dlg(n  +*(-»]</f, 

(7b) 

provided  this  limit  exists.  Applying  Eq.  7  to  average 
the  TM  cross  section  <  ID2)™  using  Eqs.  3  and  5, 
one  obtains 

<  17^10™  =  j  ( kaY(m  -  1)2(1  -  „),  (8) 


where  the  subscript  E  denotes  the  exact  result,  and  v 
■  4a/ L  is  the  packing  density.  The  term  linear  in  v 
comes  from  the  factor  <  w )  vt*  >  in  the  absolute  square 
of  Eq.  5  and  represents  mean  destructive  interference 


</U1>  = 


-  pM/4,  n  +  f  =  1 
•  VMV48,  n  +  f=  2,  (II) 
0,  n  +  t  -  4 


where  we  have  assumed  that  in  =  I  +  6  with  6  <  < 
1,  so  that  M  <  <  I .  The  terms  that  arc  second  order 
in  M  are  retained  because  the  inaccuracies  of  the 
stochastic  variational  results  first  appear  at  that 
order. 

The  exact  result  for  <  IT,  I2  >n  in  the  limits  (mka, 
ka,  and  M)  —  0  may  now  be  expressed  by  applying 


16 1.  S.  Cradshteyn  and  I.  M.  Ryzhik,  Table  of  Integrals,  Senes, 
and  Products,  Academic  Press,  New  York  (1965). 
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Eqs.  11,  9,  and  4  to  (he  average  absolute  square  of 
Eq.  6,  giving  (to  the  fourth  order  in  M) 

<  it,  1 2  >1fc  =  y  (kaYM1  [cos20,(l  -  v) 

-  Mv  cos  0,  cos  (2/3  -  dj 


As  in  the  TM  case,  the  (1  -  v )  term  describes  scatter¬ 
ing  by  two  independent  cylinders,  modified  by  inter¬ 
ference  effects  (the  -  i>  term).  The  remaining  terms  in 
Eq.  12  are  due  to  multiple  scattering. 


+  M1  ^  [2 cos2 0,  +  cos2  (2/3  -  0,)]] . 


(12) 
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4.0  VARIATIONAL  AND  BORN  APPROXIMATIONS 


Vector  stochastic  variational  expressions  were  re¬ 
cently  derived  for  general  conducting,  dielectric,  and 
magnetic  scatterers. 5,6  In  order  to  investigate  the  ex¬ 
tent  to  which  these  expressions  can  account  for  polar¬ 
ization  and  interactive  scattering  phenomena,  a 
stochastic  variational  approximation  to  the  two-Ray- 
leigh-cylinder  random  scattering  problem  addressed 
in  the  previous  sections  will  now  be  calculated.  The 
Born  approximation  for  this  problem  is  computed  as 


an  intermediate  result  that  the  stochastic  variational 
calculation  improves.  We  are  also  able  to  evaluate 
the  direct  average  of  the  deterministic  variational  ex¬ 
pression  for  this  simple  random  scattering  problem  in 
the  limit  of  nearly  transparent  Rayleigh  cylinders, 
where  m  —  I.  The  various  approximate  solutions  are 
compared  with  the  exact  solutions  in  the  Summary 
and  Conclusions. 


4.1  Variational  Formulation  for  Infinite  Cylinders 


The  variational  expressions  derived  in  Refs.  5  and 
6  are  for  scatterers  localized  in  three  dimensions  in 
which  the  scattered  field  is  a  spherical  wave  in  the 
far-field  limit.  For  the  infinite-length,  homogeneous 
cylinders  of  current  interest,  the  scattered  field  E„  is 
a  cylindrical  wave  in  this  limit  and  is  expressed  as  E„ 
—  A  T e‘k'°  /vV0  for  r„  >  >  f  (see  Fig.  1).  An  integral 
wave  equation  for  two-dimensional  scatterers  may  be 
obtained  from  the  three-dimensional  results' 6  by  rec¬ 
ognizing  that  E,  E,  and  U  =  k2(m2  -  I)  are  all  inde¬ 
pendent  of  the  z  coordinate  in  Fig.  1  for  a  normally 
incident  plane  wave.  Using  the  relationship17 


where  cs  denotes  integration  over  the  scatterer  cross 
section,  E,  is  the  incident  plane  wave  with  polariza¬ 
tion  e,  and  propagation  direction  k,,  the  two- 
dimensional  vector  f  =  xx  +  yy,  and  the  two- 
dimensional  dyadic  Green’s  function  is 


GJs'(f,n  =  Py\l+  vv/Ar:] 

x  -Hi"ikR,)  -  us) 

4 


G„( r,  r  t 


'H„  *»(k/?2) 


(13) 


with  k2  =  k2  +  h2  and 

R2  =  sA  (*-*') 2  +  ( y-y ')2, 

one  can  integrate  the  three-dimensional  results'6 
over  tne  z  coordinate  to  obtain 


E(f)  =  E,  + 


(  c<2un 


(UE<njdS', 


(14) 


Flere  /  is  the  three-dimensional  unit  dyadic,  H{‘" 
(kR2)  is  the  0-th  order  Hankel  function  of  the  first 
kind,  PV  indicates  principal  value  integration  when 
the  first  term  in  Eq.  15  is  in  the  integrand,  and  /.  is 
the  depolarization  dyadic  that  depends  on  the  shape 
of  the  infinitesimal  excluded  region  used  in  com¬ 
puting  the  PV  integral.  (We  note  that  the  total  in¬ 
tegral  involving  G„2’  is  independent  of  this  shape."1) 
In  the  far  field,  because  f  and  f '  are  never  coinci¬ 
dent,  Eq.  14  can  be  used  to  express  the  scattering 
amplitude,  T  =  e,  •  T,  as 


T  = 


(16) 


|7L.  B.  Felsen  and  N.  Marcuvitz,  Radiation  and  Scattering  of 
Waves,  Premice  Hall,  Inc.,  Englewood  Cliffs,  New  Jersey,  p. 
635(1973). 


i8A.  D.  Yaghjian,  “Electric  Dyadic  Green’s  Functions  in  the 
Source  Region,"  Proc.  IEEE  (A,  248-263  (1980). 
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with 


where 


N,  =  dS'is  •  [/-*,*,]  •  [UE(n]e  *‘  r  , 

(17) 

where  e,  specifies  the  scattered  wave  polarization  of 
interest,  ks  is  the  unit  vector  in  the  scattered  direction 
that  makes  an  angle  0a  relative  to  the  incident  direc¬ 
tion  in  Fig.  1,  and  integration  is  over  the  cross-sec¬ 
tional  area  of  the  scatterer. 

The  procedures  outlined  in  Refs.  S  and  6  may  be 
applied  to  Eq.  14  to  derive  a  stochastic  variational 
principle  for  T.  As  for  the  three-dimensional  case, 
this  results  in  an  expression  of  the  form 


N^.k.^X)  =  (19) 

and 

Z>=[  dS'  UE(f')  ■  E(f')  -  (  ds(  dS'  I  t/E(f)  ] 

*  cs  J  cs  J  cs 

■  G0<*(f,n  •  [(/E(f')l  •  (20) 


<  17V  I2  > 


J_  <)A,  l2><! JV2I2> 
16  <  IZ?  1 2  > 


(18) 


4.2  Evaluation  of  the  Integrals  and  D 


We  encountered  calculational  difficulties  in  our  at¬ 
tempts  to  evaluate  the  double  surface  integral  in  Eq. 
20  for  two  dielectric  cylinders.  Consequently,  we  em¬ 
ployed19  methods  analogous  to  those  used  by  Yagh- 
jian"  to  obtain  an  alternative  expression  in  which  D 
is  reexpressed  as 


D  *  D,  -  D,  +  D,  ,  (21) 


where 


D, 


dS'Ut(i')  ■  [(1  +  CZ/A^EO-')!  , 

(22) 


xd/Etf'WWOWi).  (23) 


l9J.  A.  Krill,  R.  H.  Andreo,  and  R.  A.  Farrell,  “A  Computational 
Alternative  for  Variational  Expressions  that  involve  Dyadic 
Green  Functions,”  IEEE  Trans.  Antennas  Propag.  AP-M, 
1003-1005(1982). 


and 


dSU E(0  -  §vdtn'  •  (cE(f') 


x  [  V  'H"\kR2)  x  /] 


+  IV'  x  imn\  X  />/«"(A:/?2) j,  (24) 


where  the  contour  integral  in  D,  is  over  the  boundary 
T  of  the  scatterer  cross  section  (cs),  and  n  is  the  out¬ 
ward  normal  at  the  boundary.  Although  Eq‘.  21  is  not 
as  compact  as  Eq.  20,  its  evaluation  is  straightfor¬ 
ward  because  ail  of  the  singularities  are  integrable. 

For  Rayleigh  dielectric  cylinders  in  which  ka  <  <  1 
and  mka  <  <  1,  it  is  appropriate  to  approximate  the 
total  original  and  adjoint  fields  appearing  in  the  in¬ 
tegrals  of  the  variational  expression  by  the  cor¬ 
responding  fields  that  would  occur  in  the  absence  of 
the  cylinders,  i.e.,  by  the  original  and  adjoint  plane 
waves5  4  E  *  Ae,  exp  (tk,  •  f)  and  E  *  Aes  exp  ( -  <k, 
•  0,  respectively,  where  f  m  r  -  zi.  Using  these  trial 
approximations,  one  finds 


/V,  =  Nt  =  A(is  ■  -  1)  £  f  e”  . 

j  icsi 

(25) 
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where  the  sum  is  over  cylinders  j  =  ±  1  and  the  in¬ 
tegral  is  over  the  cross  section  of  the  y'-th  cylinder. 
Further,  we  have  used  the  definitions  a  ■  It,  -  k, 
and  U  =  k2(m2  -  1)  for  f'tcsJt  and  U  =  0  for 
t'feSj  .The  integral  is  evaluated  by  first  expressing  the 
coordinates  for  points  interior  to  the  cylinders  as  £ 
=  f01  +  r,  where  =  [xf01  cos  0  +  y^'  sin  £]  is 
the  position  of  the  y'-th  cylinder  axis  (see  Fig.  1)  and  r 
is  the  field  point  within  the  y'-th  cylinder  relative  to  its 
axis.  Applying  this  transformation  to  Eq.  25  results 
in 


Nt  =  A(is  e,)k2(m2  -  DiriHS,  ,  (26a) 


and 


Dy  =-  A2k2(m 2  -  1)!  £  £  (  dSe  *>  f<e, 

4  ,  j  Jcsj 

■  §  dr  A'  ■  \e,  e*'  0  X  |v'W(Ar/f2)  *  f| 


+  [v'  x  (e‘  eA'  0 )]  X  //'"(*/?,)/}.  (29) 


where 

S,  =  Ee'°  f0>  =  e'"  f"’  (1  +  e‘“  f), 

J  (26b) 


where  Ty  denotes  the  circumference  of  the  y-th 
cylinder.  By  applying  the  coordinate  transformation 
f/  =  VJ>  +  r  and  evaluating  the  double  curl  in  Eq. 
28,  D2  reduces  to 

D,  =  -A2k2(m2  -  I)*V5  •  «,) 


a  *  lor  I  , 


and  we  have  used  the  fact  that 


E 


'dS  =  xa1 


r  2y,(aa)  i 

I  M  J 


J  a— 0 


‘X«r 


Applying  these  trial  approximations  to  the  denom¬ 
inator  integrals,  Eqs.  22  through  24,  results  in 


D,  =  A2(e,  •  £,)k2  (m2  -  \)m2 


x  £  e“  (jdS'  =  m2NiA  ,  (27) 

i 

D2=-A  A2k2(m2  -  l)2  E  E  i  dS 

4  ,  , 


ds\  dS'  e  *  e**  * 

CS,  ^  CSj 


Hl‘\kR2) 

(30) 


since  e,  -k,  =0.  By  applying  the  addition  theorem 
(e.g.,  see  Watson20)  for  the  cases  i  =  j  and  i  *  j  in 
Eq.  30  and  evaluating  the  resulting  integrals  in  the 
Rayleigh  limit,  Eq.  30  is  found  to  be  proportional  to 
(ka)* .  This  should  be  contrasted  with  the  other  in¬ 
tegrals  N,,N2,Dt,  and  (as  will  be  shown)  D} ,  which 
are  of  the  order  (*a)2 .  Thus,  D2  does  not  contribute 
to  the  variational  calculation  in  the  Rayleigh  limit. 

Evaluation  of  the  D ,  integral,  Eq.  29,  remains. 
This  integral  arises  from  the  discontinuity  in  refrac¬ 
tive  index  at  the  surface  of  the  scatterer  and  vanishes 
in  the  case  of  soft  scatterers.  **  (Of  course,  D2  would 
be  more  complicated  for  soft  scatterers.)  We  will 
show  that,  for  hard  scatterers,  £>,  accounts  for  all  the 
multiple  scattering  and  contributes  to  the  polariza¬ 
tion  dependence  in  the  variational  approximation. 
We  rewrite  Eq.  29  in  a  form  more  convenient  for  cal¬ 
culation  by  using  a  vector  identity  and  by  translating 
the  cylindrical  coordinate  systems.  This  gives 


•  [v'  x  V'  x  (e,  e*'  )]  x  H"'(kR2), 

(28) 


“G.  N.  Watson,  4  Treatise  on  the  Theory  of  Bessel  Functions, 
2nd  ed.,  Cambridge  University  Press,  New  York,  Chap.  XI 
(1962). 
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Dy-'-A'k r2(/n2  -  1)’ £  £e-*WlV, 

i  J 

n*v<fe 'rf*'e 
0 

•  [f  •  [vflTOW,)]* 


-  f(e,  •  VH‘"{kR2))  -  Hf  ±  k,)lH<"(kR2) 


+  Hf  -  e‘)k,HA"(kX2)}\^a.  (31) 


where  p  is  the  radial  cylindrical  unit  vector.  When  j 
=  j>  H2  =  [a2  +  t'2  -  2ai'  cos  -  <t>')]  Vl,  and 
when;  *  j.  Rz  =  [P2  +  a2  -  2Pa  cos  (0  -  p)  ] 7 , 
whereP  =  [f2  +  t'2  -  2ft'  cos  (x  -  <*>'  -  0))v- . 
These  geometrical  parameters  are  illustrated  in  Fig. 
2.  The  addition  theorem  is  applied  once  for  i  =  j  and 
twice  for  ;  *  j,  and  the  gradients  on  Ha(kR2)  are 
then  taken.  The  resulting  series  of  integrals  is 
straightforward  to  evaluate,  and  upon  retaining 
terms  that  are  the  lowest  order  in  ka,  we  find 


Dj  =  i42ir(*a)2(m2  -  1) 


x  (  («2  -  l)[(y  •  £,)(y  -  it)/2  -  i,  •  $]  S, 


-  j  (w2  ~  ma)2ff!"(kM  ■£,)£, 


■  [x  sin  (20)  -  y  cos  (20)  ]  S2  ] ,  (32a) 


Figure  2  —  Illustration  of  parameters  that  relate  to  fi2  us¬ 
ing  the  cosine  law.  When  /  =  /,  t  and  are  within  the 
same  cylinder;  however,  when  i  *  j,  t  and  i  lie  in  different 
cylinders. 


where 


w  £  e*>  r01  e*.  r°  ’ 


=  e"  f'"(e  *<  (  +  e*>  f),  (32b) 


and  S,  is  defined  in  Eq.  26b.  Observe  that  (kaY 
Hz"(kfi  in  Eq.  32  is  of  the  same  order  as  c^p  in  the 
exact  TE  solution  Eqs.  4  and  6.  For  f  —  2a,  this  coef¬ 
ficient  is  of  the  order  (ka)1 . 


4.3  Evaluation  of  the  Averages 

The  integrals  IN,  I2,  IN2 12,  and  ID,  I2  depend  on  Using  Eq.  7  to  average  these  integrals  over  random 
the  cylinder  separation  f  through  the  phase  term  (but  nonoverlapping)  separations,  we  find 


IS,  I2  =  2  +  2  Real  (e*  f) .  (33)  <  IN,  l2>  =  IA  l22ir 2(*e)4(m2  -  1)2(^  •  ^)2(l  -  r) 

(34) 
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and,  from  Eqs.  25  and  27, 


to  first  order  in  ».  These  results  lead  to 


<IN,I2>  =  <IMI2>  =  <ID,l2>/(m4/(2).  (35) 

The  final  step  needed  to  obtain  <  IT,  I  >-’  from  Eq. 
18  is  to  complete  the  evaluation  of  (  IDI2).  From 
Eq.  21  and  the  fact  that  D,  does  not  contribute  in  the 
Rayleigh  limit, 

<  l£>l2  >  =  <  \D,  I2  >  +  2  Real  >  +  <!£>,  I2  > . 

(36) 


<  !/>!->  =  A*r{nr  -  l)2(*a)4  [2(1  -  p)[(m  -  1): 


*  ( y  ■ «,)(«,  • «,)  +  («,  • «,)’  ] 


The  averages  appearing  in  the  latter  two  terms  of  this 
equation  can  be  expressed  in  the  general  form 


+  ^  (m-  -  l)2f(y  <, )  «,  •  [  Jr  sin  (2d)  -  _ycos(2d)]] 


</({)>  =  Hka)*‘2>  4 (Hr f)//,,2>  (ft  ft  cos  (>  ft) 


-  p(>  *,)«.  •  I J?  sin  (2d)  -  y  cos  (2d)  1 


2<fta)42''  » 
L 


x  ('  //J'""(*f)//J2*  (*f)  cosyfrff, 

J  2* 

(37) 

as  can  be  shown  by  direct  substitution  of  Eqs.  31  and 
27  (with  Eq.  26)  into  Eq.  36.  In  particular,  Eq.  37 
must  be  evaluated  for  the  cases  n  +  t  «  (0,1,2)  and  y 
=  [0,  *  cos  d,  ft  cos  (8,  -  d),or  ft(cosd  ±  cos(f  - 
d»|.  In  the  Rayleigh  limit,  the  evaluation  of  Eq.  37 
can  be  accomplished  using  procedures  analogous  to 
those  used  to  obtain  the  average  of  Eq.  9,*  with  the 
result 


(far)4  n  +  f  =  0,  y  =  0 

-r(ftu)4  n  +  (  =  0,  y  *  0 

</(()>  =  {  (-IKOr/ir) {kaV  n  +  (  =  I 
(/)"  '(r^XftaV  i  +  f=2 

(38) 


*An  alternate  evaluation  procedurL  1  pital's  rule,  in 

given  in  Ref.  21.  This  alternative  c.  i-  ■  J  to  evaluate  the 
average  of  Eq.  9. 

21 E.  P.  Cray,  R.  W.  Hart,  and  R.  A.  Farrell,  A  Variational  Ap¬ 
proximation  for  the  Scattering  of  Scalar  Waves  by  Stochastic 
Surfaces,  JHU/APL  TO  1322(1979). 


I-  1,-“'  •  '  ♦  Wr-I>«  •<)])  • 

(TO 

The  final  step  is  to  insen  Eqs.  34.  35,  and  39  into 
the  variational  expression  Eq.  18.  Expanding  the  re¬ 
sult  in  powers  of  v  and  retaining  the  first-order  term 
yields 

<ir,l2>,M  =  y(*«)4(m'  -  l)-'(l  -r)  (40) 

and 

<  I T,  I2  >"  =  y  (fta)4*#2  [cos2  9.(1  -  e) 

-  Mr  cos  #,  cos  (2d  -  0, ) 

-  A#2~  cos’  (2d  -  9.)},  (41) 


where  M  ■  (nr  -  I  )/(/*•’  +  I);  and  we  have  used 
the  facts  that 


{Z  for  TM  waves 

(♦2a) 

y  for  TE  waves 
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and 


where 


.  Z  for  TM  waves 

es  = 

— jrsinO,  +  y  cos  0,  for  TE  waves 

(42b) 

This  completes  the  evaluation  of  our  vector 
stochastic  variational  principle.  Next,  we  obtain  the 
noninvariant  Born  approximation  aH  the  direct 
average  of  the  Levine-Schwinger  dete  unistic  varia¬ 
tional  result. 

The  noninvariant  Born  approximation  <  \Tb  I2  >  is 
obtained  by  using  the  plane  wave  trial  function  in  the 
expression* 


It  should  be  noted  that  D  does  not  depend  on  M  and 
is  of  order  (Urn)2  except  when  the  cylinders  are  nearly 
touching  where  O  is  of  order  unity.  Using  Eq.  47,  we 
can  write  the  direct  average  of  Eq.  46  as 


<  irH  I-  >  *  <  I  A',  I-  >/(!6  IA  I2) .  (43) 

From  Eqs.  34  and  42,  the  TE  and  TM  solutions  are 

<IT,I2>™  =  ~{ka)4fm2  -  l)-(l  -  v) 

(44) 


(  IN, /A  I2  >/ 16 
J  [<  I N,/(/l(m2  +  l))l2> 

-M<r,>  +  Kf  <r:>] 


for  theTM  case 
(49) 

for  the  TE  case 


and 

<IT„IJ>TE  =  ^(kaYim2  -  l)2  cos2  9,(1  -  v) . 

(45) 

The  standard  Levine-Schwinger  form  of  the  varia¬ 
tional  principle  is 


I7"ts'2 


l/V,  I2  IN,  I2 

16  ID  1 2 


(46) 


For  nearly  transparent  Rayleigh  cylinders,  the 
denominator  integral,  D,  can  be  expanded  in  powers 
of  M.  In  particular,  recognizing  that  D,  =  0  in  this 
limit  and  combining  Eqs.  27  and  32, 


D  =  D,  +  D, 


where  IN,  /(A  {in'  +  1)) 1 2 ,  r, ,  and  I\  are  each  pro¬ 
portional  to  A#2,  and  terms  of  order  Kf  and  higher 
have  been  neglected.  In  particular. 


T,  =  i~  {ka)*M:  cos 0,  cos(2d  -  0. ) 

4 

x  [«,"'<* osrs,  -  «,'2,<*f)s,sn  (50) 


and 

r,  =  K4\f  (y  )  COS:  (2d  -  0) 

x  [|  (*a) HY  <*f)S,  f 


N,  A  for  the  TM  case 

N,A  (  —  ) 1 1  +  MO)  fortheTEcase 

2  (47) 


-  <*(i)j («("(*j a ):  srs//s, 

-  (ka),(f/P(kny  S,S,*2/Sfj.  (51) 
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We  note  that  S*S2  and  S2/S,  depend  on  the  cylin¬ 
ders'  positions  only  through  the  axis  separation  f. 
Therefore,  the  averages  appearing  in  Eq.  49  can  be 
evaluated  using  Eqs.  37  and  38.  The  result  is 


r 


T2 

8 


(ka)\m2-  1)2(1  -  y) 


for  the  TM  case 


<I7VSIJ> 


{ka)*Ml  £cos2  61,(1  -  *)  (52) 

-  My  cos  6,  cos  (2d  -  0, ) 

+  Af2  -  cos2  (2d  -  0,)]  J  for  theTE  case 
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5.0  SUMMARY  AND  CONCLUSIONS 


The  recently  developed  vector  stochastic  varia¬ 
tional  principle  for  electromagnetic  wave  scattering 
has  been  reexpressed  in  terms  of  the  scalar  Helm- 
holtz-Green’s  function  and  evaluated  for  a  test  case 
of  two  randomly  separated  dielectric  cylinders.  This 
study  examines,  for  the  first  time,  the  extent  to  which 
this  principle  can  account  for  polarization  and  inter¬ 
active  phenomena  even  though  the  approximation 
that  it  improves  does  not.  This  has  been  accom¬ 
plished  by  calculating  the  exact  solution  and  compar¬ 
ing  the  result  with  the  variational  approximation. 
The  two  randomly  separated  cylinders  have  parallel 
axes  and  are  infinitely  long.  The  cylinders  are  not 
allowed  to  overlap,  i.e.,  they  are  hard,  and  their  di¬ 
ameter  is  small  relative  to  a  wavelength.  The  varia¬ 
tional  and  exact  results  may  also  be  compared  with 
the  Born  approximation,  which  is  obtained  as  an  in¬ 
termediate  step  in  calculating  the  variational  result 
when  simple  plane  wave  trial  functions  are  used. 

Comparison  of  Eqs.  8,  40,  and  44  indicates  that 
both  the  Born  and  the  stochastic  variational  approx¬ 
imations  give  the  exact  result  for  the  average  TM  dif¬ 
ferential  scattering  cross  section.  As  the  exact  TM  so¬ 
lution  contains  no  multiple  scattering,  this  conclu¬ 
sion  is  not  surprising. 

Comparison  of  Eqs.  12  and  45  indicates  significant 
discrepancies  between  the  Born  approximation  and 
the  exact  result  for  the  mean  TE  cross  section.  First, 
the  Born  approximation  accounts  for  the  angular  de¬ 
pendence  of  the  TE  solution  but  does  not  give  the 
correct  refractive  index  dependence  in  its  amplitude, 
i.e.,  (m2  -  l)2/8  for  the  Born  versus  \f  /2  s  [(m2 
-  l)/(m2  +  1 )  ] 2  /2  for  the  exact  solution.  This  dis¬ 
crepancy  arises  because  the  Born  approximation  neg¬ 


lects  the  geometric  polarizability  (resulting  from  the 
cylindrical  shape)  in  approximating  the  fields  inside 
the  cylinders.  Second,  as  expected,  the  Born  approx¬ 
imation  includes  interference  but  does  not  account 
for  TE  multiple  scattering. 

The  stochastic  variational  improvement  (Eq.  41 )  of 
the  Born  approximation  for  this  TE  scattering  cor¬ 
rects  for  the  geometric  polarizability.  In  addition,  it 
contains  a  multiple  scattering  contribution  that  is 
correct  through  the  lowest -order  term  in  the  variable 
A/,  i.e.,  through  A/’.  Thus  the  stochastic  variational 
correction  factor  1/4  I2  <  IN,  I2  >/<  ID  1 2  >  significant¬ 
ly  improves  the  Born  approximation  to  give  the  po¬ 
larization  and  multiple  scattering  dependences  ac¬ 
curately  (even  though  these  effects  are  missing  from 
the  simple  Born  trial  field). 

The  direct  average  of  the  Levine-Schwinger  deter¬ 
ministic  variational  expression  was  obtained  in  the 
limit  of  nearly  transparent  Rayleigh  cylinders.  Thai 
approximation  also  agrees  with  the  exact  result 
through  the  terms  of  order  A/' . 

As  with  the  stochastic  variational  principle,  the  M‘ 
term  is  incorrect.  (This  discrepancy  was  uncovered 
subsequent  to  the  publication  of  Ref.  22.)  Neither 
variational  result  reproduces  the  cos2  0,  dependence 
found  in  the  M*  term  of  the  exact  result,  and  both 
give  an  incorrect  value  for  the  coefficient  of  the 
cos2(20  -  8. )  term. 


22J.  A.  Krill,  R.  H  Andreo,  and  R.  A  Farrell,  "Variational  C'alcu- 
laiions  of  Electromagnetic  Scattering  from  Two  Randomly 
Separated  Rayleigh  Dielectric  Cylinders."  J.  Opi.  Soc.  Am.  73. 
408-410(1983). 
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